Preheating with Non-Standard Kinetic Term 
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We consider reheating after inflation in theories with a non-standard kinetic term. We show that 
the equation that exhibits parametric resonance is the Hill-Whittaker equation, which is a particular 
case of the more general Hill equation. This equation is in general more unstable than the Mathieu 
equation, such that narrow resonance preheating may be efficient in theories with a non-canonical 
kinetic term. 

PACS numbers: 98.80.Cq 



I. INTRODUCTION 

Reheating is the period after inflation during which the 
energy stored in the inflaton is transferred to the matter 
and radiation of which the universe is made up. During 
this stage, the inflaton oscillates about its minimum and 
decays into relativistic particles. Through coupling with 
other fields, these oscillations can give rise to parametric 
instability and lead to an explosive growth inparticle 
density. This phase, called preheating 0, 0, y, d, ||, 
plays a crucial role in the dynamics of reheating of an 
inflationary universe. 

It has been shown 0, 0] (see also [1, that para- 
metric amplification of entropy modes during reheating 
can lead to a change in the curvature power spectrum on 
scales larger than the Hubble radius (but smaller than 
the horizon) without violating causality. Indeed entropy 
perturbations act as a source for large scale curvature 
perturbations. The effect is negligible for single-field in- 
flation - where entropy perturbations are suppressed on 
super-Hubble scales - but can be large in the case of mul- 
tiple field inflation [ic| . The study of this mechanism in 
multiple field models demands attention since it can sig- 
nificantly modify the usual inflationary picture. 

It has been shown [ll[ that for the simplest preheating 
model, the contribution to the large scale power spec- 
trum of curvature perturbations from entropy modes is 
negligible. However, very little work has been done to 
this date to show that this is also the case for more com- 
plicated scenarios. With the large number of inflationary 
models proposed today, one can expect the study of this 
mechanism to become a powerful tool. 

In this paper, we study preheating in two-field infla- 
tion models {</>,%} with a non-standard kinetic term for 
X- Such a kinetic term arises inevitably in string-inspired 
models where the Kahler potential K. has a non-trivial 
geometry. Examples of this are brane inflation models, 
where the inflaton is the distance between two branes 
in higher dimensional spacetime, and modular inflation, 
where inflation is caused by the dynamics of the defor- 
mations of a six dimensional Calabi-Yau manifold (see 
e -g- EH EH EH El| for reviews on string inflation 
models). 

Using the general form of a two-field Lagrangian with 
a non-standard kinetic term introduced in [l7| (see also 



(18j|). we show that the entropy fluctuations during re- 
heating obey the Hill-Whittaker equation, a modified ver- 
sion of the Mathieu equation (the equation that exhibits 
parametric instability). Indeed, the non-canonical part 
of the kinetic term enters as an oscillatory damping term 
into the equation, modifying its stability behaviour. We 
show that the Floquet theorem applies for this new type 
of equation such that the modes Xk undergo paramet- 
ric resonance for a wide range of parameters. Regions 
of instability appear as bands in the parameter space. 
We show that they are larger in the case of the modified 
Mathieu equation. Therefore preheating is more likely to 
happen in models with non-standard kinetic terms. 



II. PREHEATING 

Consider a theory of two interacting scalar fields with 
non-canonical kinetic term [13, EH 



S = Jd 4 x^g~[^fR 



\{d^){d^) 



-(d^x)(d^x)-V(^ X )}, 



(1) 



where Mp is the reduced Planck mass and b((f>) is a non- 
trivial function of <f> that renders the kinetic term non- 
canonical. This typ e of action appears in the Roulette 
inflation model I l 9 f . It also appears in generalized Ein- 
stein theories |20l.l2ll] . In a Friedmann-Robertson- Walker 
(FRW) background 



-dt l 



a(t) 2 dx 2 , 



(2) 



with a(t) denoting the scale factor, the equations of mo- 
tion are 



3H<t> + V (j> = 



X+(3H + 2b^) X + e~ 2b V x = 



H 2 = 



1 



3Mp 



1 



-,2b 



H = - 



2Mp 



e 2b X 2 



(3) 
(4) 

(5) 
(6) 



Note that for 6 = the equations of motion for the 
fields 4> and x reduce to the Klein-Gordon equations and 
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Eqs.©-© reduce to the well known Friedman equations 
in the presence of scalar field matter. 

Consider a model of chaotic inflation with a symmetry 
breaking potential 



f 



1 



1 



V(<f>, x) = 2™ W> - CT ) + 2 m x* + W (7) 

where is the inflaton and x is the reheat field. We 
start with a vanishing x background field and will study 
the growth of the perturbations $x which begin as initial 
quantum vacuum fluctuations. Inflation occurs for \<f> — 
crj > Mp. In the case of a standard kinetic term (b — 
0) and in the limit m 3> m x , this model gives efficient 
preheating p], Q . We shall investigate the effect on this 
system of a non- vanishing &(</>). 

After the shift <fi — a — > <fi, the potential becomes 



v(<f>,x) 



2i2 i 1 2 2 , 1 2j.2 2 2i 2 

m <p + -m x x + -g <p X + g WX 



,222 

+ 2 9 ° X ■ 



(8) 



When <fi ~ Mp, inflation ends and the inflaton starts 
oscillating about the minimum of its potential 



(j)(t) = $(t)cos(mt) , 



(9) 



where $(f) is the amplitude of oscillation. This ampli- 
tude decreases as ~ 1/t due to Hubble friction. In the 
case when the period of oscillations is small compared to 
the Hubble time, i.e. m ^> H, we can neglect the Hub- 
ble friction as a first approximation: preheating happens 
on a time scale much shorter than the Hubble scale and 
we expect other effects - such as back-reaction from the 
newly created x particles - to come into play long before 
Hubble friction is felt. 

Consider small amplitude oscillations $ <C u, such 
that the term \g 2 (f) 2 x 2 m the potential can be neglected. 
Then, the equation for perturbations 5\ is 

• • k 2 

5x+2b5x+—+e~ 2b [g 2 a 2 + m\ + 2g 2 a<i>cos(mt)] Sx = 

(10) 



To first order in (— ), the oscillatory friction term is 



4> 



msin(mt) 



(11) 



and the prefactor of the last term on the left hand side 
of lUDt takes the form 



-26 



where b° {rj>/a) 



-2b 



i 1 ■- 2b (<t>M (?)"-("'/) 



(12) 



db 



and b n 



6=0 



6=0 



Performing the change of variables mt 



2z, and ne- 



glecting 0((— Y) contributions, we get 
S X " + 2p ■ sin(2z)5x' + [A k - 2qcos{2z)]5x = , (13) 



where a prime denotes the derivative with respect to z 
and 



A 4fc 2 , 4e" 26 » ( ., , ,^ 



(14) 



Eq. (|10p is reminiscent of the well known Mathieu equa- 
tion, but with an additional oscillatory friction/anti- 
friction term with coefficient 2p. Defining the new func- 



tion f(z) 



-pcos(2z)/2 



Sx, Eq.([T3"j) gives 



f"+ A fe -^--2(g + p)cos(2z) + ^-cos(4z) / = 0, 

(15) 

which is known as the Whittaker-Hill equation and is a 
particular case of Hill's equation [22|, [23|. The Floquet 
theorem applies for this type of equation and implies that 
any solution can be written in the form 



/(*) = ae" 2 0(z) + /3e-» z (f>(-z), 



(16) 



where /i e C is called the Floquet exponent, <fi(z) is pe- 
riodic in z with period 7T, and a, (3 are real constants. If 
\x is imaginary, the solution is stable whereas jU real gives 
exponential instability. Note that the Floquet exponent 
fi - which determines the rate of exponential growth of 
solutions - is the same for Sx(z) and f(z). 



III. STABILITY OF THE WHITTAKER-HILL 
EQUATION 

The Floquet theorem implies that the Whittaker-Hill 
equation can exhibit exponential instability. A method 
for calculating the resonance was derived by Whittaker 
in [24( for the full Hill equation, which includes the 
Whittaker-Hill equation as a particular case. For clar- 
ity, we present the details of this method. We start by 
writing a solution of the form 



y(z) = e^cf>(z) = e^X> r 



_2riz 



(fi-\-2ri)z 



(17) 



and plug it into Eq. (fH)|) . We obtain the following recur- 
sion relation for the coefficients: 



(c 2 (r+2) + C 2(r _2)) TT 



p 2 /A 



(ifi — 2r) 2 — a 



+ (c 2 (r-l) + C 2 (r+1)) TJ 



q + p 



{ill — 2r) 2 — a 



+ c 2r = (18) 



3 



where we use the abbreviation a for A^, which can be 
written as 

c 2(r-2)£2r + C 2 ( r -l)l2r + C 2r + C 2 ( r+ l)72r + C2(r+2)£,2r = 0, 

(19) 

where 

6, S „ p2 Ji • (20) 



72r 



(z/i — 2r) 2 — a 



(j/i — 2r) 2 — a 



The value of ifi which solves this recursion relation is that 
for which the determinant of the coefficients vanishes 



AM 



£-2 7-2 1 7-2 £-2 

£o 7o 1 7o fo 

6 72 1 72 6 



0. 



Let us investigate the singularities of this determinant in 
the complex plane. For y/a not an integer, the determi- 
nant has simple poles at 



ifi = 2r ± y/a. 



(21) 



There are pairs of poles at distance y/a of every point 21 
on the real line (I a real integer). Let us show that the 
residues at these poles alternate in sign and are periodic. 
This will allow us to construct an entire function Q(z) 
which will help us solve the recursion relation. 

Theorem 1. 

(i) A(ifj) is periodic in the reals with period 2m, 

where m is a real integer, 

(ii) A(ifj) is even in the reals about every point 2m. 

(22) 

Proof. One can check that the shift ifj, — > ifj, + 2m, (m a 
real integer) induces the change 



72r — * 72(r-m) 
?2r — * ?2(r-m) 



(23) 



such that the rows of the determinant are all shifted up 
by an equal amount. This does not change the value of 
the determinant of the infinite matrix. To prove (ii), one 
can check that the change ifi — > —ifj, induces the shift 



72r -> 7-2r 
£,2r — > £-2r 



(24) 



such that the rows of the determinant are nipped about 
the row with the 7o's, i.e. 



AM 



6 72 1 72 6 

Co 7o 1 7o Co 

£-2 7-2 1 7-2 £-2 



0. 



From the method of cofactors one can check that this is 
equal to A(ifi). □ 



Corollary 2. For any even integers n f2, i/ie residues 
of A(i/i) satisfy 

(i) Res(A(ifj,),n + \fa) = Res(A(ifi), r 2 + y/a), 

(ii) Res(A(ifi) : n — \/a) — Res(A(ifi),r 2 - y/a), 

(Hi) Res(A(ifi), r% + y/a) = —Res(A(ifi), n — yfa), 

(25) 

where Res(A(ifi), z) stands for the residue of A(ifi) at 
point z. 

Proof. The proof follows directly from Theorem [TJ Let 
us perform a Laurent expansion about each of the poles. 
In some neighbourhood around i\x — y/a, we have 



AM = 



(ifi - y/a) 



+ C + d(ifi - y/a) + 



(26) 



whereas, in some neighborhood around ifi = —y/a, we 
have 



AM 



D-i 



(ifi + y/a) 



+ D + D x (ifi + y/a) + 



(27) 



We see that A(ifi) is even if and only if C-\ = — -D-i, 
Co = Dq, C\ = —D\... The constants C_i and D-\ are 
the residues at the poles. □ 

Let us introduce a function A which has the same an- 
alyticity structure as our determinant. 



AM 



1 



cos(irifi) — cos(TTy/a) 



(28) 



This function is periodic ifi — > ifi + 2m and has simple 
poles at ifi = 2r± y/a, just like A(ifi). From the period- 
icity of A, we conclude that the residues Res(X, 2r + y/a) 
are the same for all r and similarly for Res(X, 2r — y/a). 
Left to prove that the residues alternate in sign, i.e. 

Res(X,2r + y/a) = -Res(X,2r- y/a). (29) 

From l'Hopital's rule, the residues can be calculated ex- 
plicitly 

lim (i^j- r ^ y/a)X = ^ — —j — t=t, (30) 
i/i^2r±v/a irsin(iTy/a) 

which shows that the sign of the residue is switched be- 
tween the two types of poles. The poles of the function 
A match the poles of A(ifi) up to a constant. Thus an 
entire (analytic everywhere) function can be constructed 



C = A(ifi) — kX, 



(31) 



where K is a constant set to cancel exactly the divergences 
at the singular points. The functions A and A(ifi) have 
limiting behavior 



lim A =0, 

Im(ifi)— >oo 

lim A(i^/) =1, 



(32) 
(33) 
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since 7 r goes to as \ifj,\ — > oo (in which case A(«/i) 
approaches the identity matrix). Therefore, 



lim 



CM 



1. 



(34) 



This implies that £ is bounded. From Liouville's theo- 
rem an entire function whose absolute value is bounded 
throughout the z-plane is constant. Thus £ = 1 through- 
out the z-plane, in particular £(0) = 1. From this we 
get 



C(0) = 1 = A(0) - /sA(0), 



(35) 



such that 



A(i(i) = 1 



(36) 



(37) 



« = = (A(0) - 1)(1 - cos(n^))- 

From (|3ip. we have 

(A(0)-l)(l-coa(7rV5)) 
cos(mfi) — cos(iTy/a) 

The determinant equation is A(i/x) = 0. It is solved for 

cos(mn) = 1 + A(0)[cos(7rVa) - 1]. (38) 

The expression for /i (which is derived by Whittaker in 
[H) follows 



Next, we note the rich structure of the parameter 
space. The symmetry about the q = axis that appears 
for the usual Mathieu equation is broken when p acquires 
a non-vanishing value. Moreover, the resonance bands 
are shifted from their usual positions (Ak = 1, 4, 9, ...) as 
they become larger. 




.,t : m. 



- 1 {l + A(0)[cos(7r^)-l]} 



(39) 



This expression has the advantage that it can be eval- 
uated efficiently using numerical methods. Indeed, the 
quantity A(0) is easily approximated by taking numer- 
ically the determinant of an n x n matrix with large n 
(the center of the matrix is at the position of the 1 in the 
row containing the 70). Since the off-diagonal elements 
all scale as 1 jr as we go away from the center of the ma- 
trix, there will be a rapid convergence of the determinant 
computed in this approximation. 

Figure [1] show the real part of the Floquet exponent \x 
in the (a, q) parameter space for different values of p. The 
numerical simulations were performed using a 500 x 500 
matrix The instability bands of the p — figure match 
exactly those of the well-known Mathieu equation, which 
demonstrates the accuracy of the procedure. The first 
observation is that for bigger values of p, the instability 
bands get larger. This indicates that inflationary models 
with non-standard kinetic terms can exhibit more effi- 
cient parametric resonance. The heuristic reason for the 
more efficient resonance is the fact that the periodic vari- 
ation in the friction/anti-friction term in the equation of 
motion contains energy which can be given to the sys- 
tem. For a more mathematical argument, let us return 
to the expression (|3"9"|) for the Floquet exponent [i. If we 
increase the value of p, then the value of the determinant 
A(0) will increase. This will then lead to an increase in 
the absolute value of \x. 



\i 1.4 



0. 

04- 



H 14-, 



I 



A, 1 



FIG. 1: The Floquet exponent fi for (from top to bottom) 
p = 2, p = 1, p = 0.5 and p = 



■5 



IV. POTENTIAL WITHOUT SYMMETRY 
BREAKING 

Consider the simple potential for chaotic inflation with 
a massless \ ne ld 



We consider the situation *— 3- < H, such that the reso- 
nance is inefficient for 6 = 0. 

Consider e~ b ° of order unity or smaller. We show that 
preheating can still happen in this case. First, note that 
p dominates over q: 



t r 1 2j2 1 1 2 j.2 2 

V = -m <f> +-g (/) x ■ 



(40) 



For this type of potential, the condition for efficient pre- 
heating (with m 3> H and in the absence of non-standard 
kinetic term) is [5j 



■9 4$4 -> Tf 

q to = -5—3 > h. 



(41) 



If condition (|4Tj) is not satisfied, the modes are red-shifted 
outside the resonance band before they have time to get 
amplified and preheating does not takes place. However, 
this condition translates to 



m -1/3/^4/3 



M P g 



m . 



(42) 



Since m must be much smaller than the Planck mass, 
but 4> starts oscillating with a value close to the Planck 
mass, this condition is easy to satisfy. 

We show that by introducing non-vanishing b we can 
relax this condition. Consider small oscillations <j> <C Mp 
about the minimum. We expand 



-2b 



-2b a 



2b °WM P ) 



where b 



(<t>/M P 



db 



d(*/Mp) 

the perturbations becomes 



:>=0 



such that the equation for 



S X " - 46? 



{2e- 
Ak 2 



2b, 



■}s x = 0, 



(44) 



^x"-46^ /Mp) ( — )sin(2z)5 X '- 



4fc 2 

— = 0. (45) 



such that as a first approximation, the parameters of the 
modified Mathieu equation can be written in the form 



A k = 
9 = 

P = - 



4fc 2 



2b °<P/M F 



( — 

\M P 



(46) 
(47) 

(48) 



The graphs of Fig. Q] shows that the modified Mathieu 
equation exhibits resonant behavior on the q = axis for 
p non- vanishing, indicating that efficient resonance may 
be possible. Figure [2] shows the Floquet exponent on the 
q = axis for different value of p. 




- p-i 

- p=0.7 
p=0.5 
p=0.2 



FIG. 2: The Floquet exponent fi on the q — axis for different 
values of p. 



after dropping the higher frequency contribution propor- 
tional to cos(6z). Two cases are of interest: first, if 

b °{4>/Mp) i^Mpj ^ 1> ^ e damping term becomes negli- 
gible and cos(4z) dominates, such that the equation can 
be turned into a standard Mathieu equation after the 
change of variables 2z — ► z. In the condition q 2 m > H 

for efficient resonance, we now have q — e~ 2h ° 9 2 ^i ■ This 
case is realized in the Roulette Inflation model, which is 
discussed in a companion paper. Since to is much smaller 
than M p , the efficient resonance condition is satisfied. 



The second case is b?^, M ^ 



In this case the higher frequency contribution cos(4z) in 
the mass term can be neglected as a first approximation. 



One mechanism which can prevent resonance from be- 
ing efficient is the redshift of the modes outside the reso- 
nance band before they get sufficiently amplified. While 
in the band, the growth is 



(49) 



where 8x0 is the initial value of §x an d At is a small time 
interval. Summing over the time intervals and taking the 
continuity limit gives 



5x-5 Xo e^^ dt , 



(50) 



with [i = n(t). The condition m 3> H implies that 
H ~ constant on time scales relevant for preheating. The 
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growth (|50p can be integrated on the A k by performing 
the change of variables dt = dAk/H. It becomes 

Sx^S X oe^^ dA ". (51) 
For the growth to be effective for some mode k, we need 



2H 



n{A k )dA k 



2H 



I > 1. 



(52) 



The integral /, which is the area under the curve shown 
in Fig. [2l is an even function of p (we computed the inte- 
gral in the approximation of setting the factor A^ 1 in the 
integrand to 1 which is a good approximation since it cor- 
responds to the center of the resonance band). Plotting 
/ for p G [—1,1] shows that I(p) can be approximated by 
a parabola on this interval (see Fig. [3]). 



I(P) 



-1 -0.8 -0.6 -0.4 -0.2 



0.2 0.4 0.6 0.8 



FIG. 3: Growth factor J as a function of p. 

A least square fit gives 

I{p) ~ 0.75p 2 . 



During reheating H 



MpV6 



and condition 



(53) 
becomes 

(54) 



This condition does not contradict our first assumption 
i'Mp) ^(4>/M P ) — f° r smau amplitude oscillations. 

This implies that preheating is efficient in this case. 

Efficient parametric resonance is obtained in this 
chaotic inflation model for a broad range of parameters. 
Indeed, the introduction of a non-standard kinetic term 
that satisfies 



e-°° < 1 
b° > ( MP 



(55) 
(56) 



gives rise to preheating whatever the value of g 2 <S> 2 /m 2 . 
This simple example teaches us that reheating in theories 



with non-standard kinetic term must be approached with 
caution since parametric resonance may happen in exotic 
ways. 

Those results hold for small amplitude oscillations. No 
direct comparison with the standard large field reheat- 
ing results are possible at this point, since the expansion 
(|44| becomes inaccurate for large amplitude oscillations. 
Note that in the case we studied here, the condition q > 1 
is easily satisfied in general and therefore we should ex- 
pect to see broad resonance. The precise effect of a non- 
standard kinetic terms during the phase of broad reso- 
nance have yet to be investigated. 



V. CONCLUSIONS AND DISCUSSIONS 

In this paper, we have studied the simplest preheating 
model with a non-standard kinetic term. We observed 
that the non-standard kinetic term introduces an oscil- 
lating friction term in the Mathieu equation that mod- 
ifies the stability behavior of the perturbations in the 
reheat field %■ The equation that is obtained is the Hill- 
Whittaker equation, which in the form in which it ap- 
pears here is more unstable than the Mathieu equation. 

Thus parametric resonance is more likely to be effi- 
cient in theories with a non-standard kinetic term. In 
Section ITV1 we showed that for the simplest preheating 
model, there exists a function b(4>) that makes paramet- 
ric resonance efficient whatever the coupling g. Thus 
one should be careful when dealing with theories with 
non-standard kinetic terms in the context of reheating as 
models that do not seem to give rise to preheating at first 
sight (because they have a small value of q) can become 
resonant due to the oscillatory friction/anti- friction term 
introduced by the non-standard kinetic term in the ac- 
tion. Note, however, that for small values of b(<fi), as is 
realized in the Roulette Inflation model, the effects of the 
non-standard kinetic term are negligible. 

As discussed in a companion paper, equations like 
those studied in this paper occur in the Roulette Infla- 
tion model. There, the perturbations in the reheat field 
X correspond to entropy fluctuations that are parametri- 
cally amplified through the mechanism discussed here. 
These entropy modes seed curvature perturbations on 
large scales, altering homogeneity and scale invariance. 
In a companion paper, we were able to put constraints 
on the parameters of the Roulette Inflation model so as 
to match the results from the COBE experiment. Similar 
work was conducted in the cases of the D3/D7 brane in- 
flation model [H and the KKLMMT model In both 
cases it was found that under certain assumptions, dur- 
ing reheating entropy fluctuations can seed a secondary 
curvature mode that dominates over the primary mode 
(the purely adiabatic linear perturbation theory mode). 

Our results show that the study of parametric reso- 
nance during reheating can be used to put constraints on 
the many inflationary models proposed today. We expect 
the power of this method to be fully appreciated in the 



7 



future. a FQRNT Team Grant. We wish to thank Francis Cyr- 

Racine for valuable conversations. 

Acknowledgments 

This work is supported in part by an NSERC Discovery 
Grant, by the Canada Research Chair program, and by 



[1] J. H. Traschen and R. H. Brandenberger. Particle 

production during out-of-equilibrium phase transitions. 

Phys. Rev. D, 42(8):2491-2504, Oct 1990. 
[2] L. Kofman, A. Linde, and A. Starobinsky. Reheating 

after inflation. Phys. Rev. Lett, 73:3195-3198, 1994. 
[3] Y. Shtanov, J. H. Traschen, and R. H. Brandenberger. 

Universe reheating after inflation. Phys. Rev., D51:5438- 

5455, 1995. 

[4] P. B. Greene, L. Kofman, A. D. Linde, and A. Starobin- 
sky. Structure of resonance in preheating after inflation. 
Phys. Rev., D56:6175-6192, 1997. 

[5] L. Kofman, A. D. Linde, and A. A. Starobinsky. To- 
wards the theory of reheating after inflation. Phys. Rev., 
D56:3258-3295, 1997. 

[6] F. Finelli and Robert H. Brandenberger. Parametric am- 
plification of metric fluctuations during reheating in two 
field models. Phys. Rev., D62:083502, 2000. 

[7] Bruce A. Bassett and Fermin Viniegra. Massless metric 
preheating. Phys. Rev., D62:043507, 2000. 

[8] F. Finelli and R. H. Brandenberger. Parametric am- 
plification of gravitational fluctuations during reheating. 
Phys. Rev. Lett, 82:1362-1365, 1999. 

[9] B. A. Bassett, D. I. Kaiser, and R. Maartens. General rel- 
ativistic preheating after inflation. Phys. Lett., B455:84- 
89, 1999. 

[10] C. Gordon, D. Wands, B. A. Bassett, and R. Maartens. 
Adiabatic and entropy perturbations from inflation. 
Phys. Rev., D63:023506, 2001. 

[11] K. A. Malik, D. Wands, and C. Ungarelli. Large-scale 
curvature and entropy perturbations for multiple inter- 
acting fluids. Phys. Rev., D67:063516, 2003. 

[12] A. D. Linde. Scalar Field Fluctuations in Expanding Uni- 
verse and the New Inflationary Universe Scenario. Phys. 
Lett, B116:335, 1982. 

[13] C. P. Burgess. Lectures on Cosmic Inflation and its Po- 
tential Stringy Realizations. PoS, P2GC:008, 2006. 



[14] James M. Cline. String cosmology. 2006. 

[15] S. H. Henry Tye. Brane inflation: String theory viewed 
from the cosmos. Lect Notes Phys., 737:949-974, 2008. 

[16] Liam McAllister and Eva Silverstein. String Cosmology: 
A Review. Gen. Rel. Grav., 40:565-605, 2008. 

[17] Fabrizio Di Marco, Fabio Finelli, and Robert Branden- 
berger. Adiabatic and Isocurvature Perturbations for 
Multifield Generalized Einstein Models. Phys. Rev., 
D67:063512, 2003. 

[18] Z. Lalak, D. Langlois, S. Pokorski, and K. Turzynski. 
Curvature and isocurvature perturbations in two-field in- 
flation. JCAP, 0707:014, 2007. 

[19] J. R. Bond, L. Kofman, S. Prokushkin, and P. M. Vau- 
drevange. Roulette inflation with Kaehler moduli and 
their axions. Phys. Rev., D75: 123511, 2007. 

[20] A. Starobinsky and J. Yokoyama. Density fluctuations 
in Brans-Dicke inflation. 1994. 

[21] A. Starobinsky, S. Tsujikawa, and J. Yokoyama. Cos- 
mological perturbations from multi-field inflation in gen- 
eralized Einstein theories. Nucl. Phys., B610:383-410, 
2001. 

[22] W. Magnus and S. Winkler. Hill's equation. John Wiley 
and Sons, New York, 1966. 

[23] K. M. Urwin and F. M. Arscott. Theory of the 
Whittaker-Hill Equation. Proceedings of the Royal So- 
ciety of Edinburgh, 69:28-44, 1970. 

[24] E. T. Whittaker and G. N. Watson. A course of modern 
analysts. Cambrigde at the University Press, Cambridge, 
1962. 

[25] R. H. Brandenberger, K. Dasgupta, and A.-C. Davis. 
A Study of Structure Formation and Reheating in the 
D3/D7 Brane Inflation Model. Phys. Rev., D78:083502, 
2008. 

[26] R. H. Brandenberger, A. R. Frey, and L. C. Lorenz. En- 
tropy Fluctuations in Brane Inflation Models. 2007. 



